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Recent experiments demonstrate that a ballistic version of spin resonance, mediated by spin- 
orbit interaction, can be induced in narrow channels of a high-mobility GaAs two-dimensional 
electron gas by matching the spin precession frequency with the frequency of bouncing trajectories 
in the channel. Contrary to the typical suppression of Dyakonov-Perel' spin relaxation in confined 
geometries, the spin relaxation rate increases by orders of magnitude on resonance. Here, we present 
Monte Carlo simulations of this effect to explore the roles of varying degrees of disorder and strength 
of spin-orbit interaction. These simulations help to extract quantitative spin-orbit parameters from 
experimental measurements of ballistic spin resonance, and may guide the development of future 
spintronic devices. 



I. INTRODUCTION 

Spin-orbit interaction is the primary source of spin 
relaxation for free carriers in semiconductors.^ ^ At the 
same time, it offers the potential to control the spin orien- 
tation of those carriers without the need for conventional 
high frequency resonance techniques.'^HS Controlling car- 
rier spins using spin-orbit interaction requires the abil- 
ity to tune its effect with external parameters such as a 
magnetic field or voltages on electrostatic gates. In the 
Datta-Das spin transistor concept, for example, the spins 
of carriers in a 2D quantum well rotate in response to a 
spin-orbit interaction whose strength can be tuned by a 
gate.iasl 

Another way that electrostatic gates can tune the ef- 
fects of spin-orbit interaction in a quantum well is by 
defining the lateral confinement geometry of a spintronic 
device. Recent experiments have shown that bounc- 
ing trajectories in gate-defined channels of high mobility 
GaAs two-dimensional electron gas (2DEG), in an exter- 
nal magnetic field, lead to rapid spin relaxation through a 
process we refer to as ballistic spin resonance (BSR).^ On 
resonance the effect of spin-orbit interaction is amplified 
by matching the bouncing frequency to the Larmor pre- 
cession frequency, and the bouncing frequency depends 
on the gate-defined channel width. Although the mech- 
anism of BSR is straightforward, it is not obvious that 
the effect should be visible for realistic parameters in a 
practical device. 

In this report, semi-classical Monte Carlo simulations 
of spin dynamics are used to test the resilience of BSR 
over a wide range of device parameters. The simulation 
models varying degrees and types of disorder, confine- 
ment potential from the electrostatic gates, and lack of 
perfect specularity on scattering off the channel walls. 
We restrict our attention to electron-doped GaAs 2DEGs 
at low temperature, where the Dyakonov-Perel' mecha- 
nism has been shown to be the dominant source of spin 
relaxation.''^ ^ A range of spin-orbit interaction strengths 
are explored in the simulation, including linear Rashba 
and Dresselhaus terms as well as the cubic Dresselhaus 



term. BSR is found to be robust over a wide range of 
experimentally-accessible parameters, and not to depend 
sensitively on specific model of disorder. 

II. DYAKONOV-PEREL' MECHANISM 

At the most general level, spin-orbit interaction cou- 
ples an electron's spin degree of freedom to its momen- 
tum. Spin-orbit interaction in III-V semiconductor quan- 
tum wells is characterized by Rashba (a) and Dressel- 
haus (/3 and 7) terms, due respectively to structural and 
bulk crystal inversion asymmetry.l^ Including both types 
of spin-orbit interaction, the spin-orbit Hamiltonian is: 

Hso = "(fcoioo'ioo - ^loofQio) + /^Ifcioofioo — ^oiofoio) 
+ 7(^010^100^100 - ^^OO^OlOCToio) (1) 

where fcioo is the component of the Fermi wavevector 
along the [100] crystal axis, aioo is the Pauli spin operator 
along the [100] axis. We note that the sign convention 
adopted for this Hamiltonian gives j3 ss —7 < k^^i >.'^ 
The linear-in-fc terms in the Hamiltonian become sim- 
pler when described along [110] and [110] axes. For the 
rest of this paper, the [110] axis is referred to as the x- 
axis; the [110] axis is referred to as the y-axis. Hgo can be 
interpreted in terms of a momentum-dependent effective 
Zeeman field, Bso, which takes the following form when 
expressed along x- and y-axes: 

Bso = — ((a - (i)kyx ~{a + P)Ky) + 0{k^), (2) 

This effective field corresponds to a spin-orbit precession 
time Tgo = fi/ gfJ-sBsoi where g = 0.44 is the Lande g- 
factor in GaAs. 

Spin precession according to Hso is coherent from a 
microscopic point of view. At a practical level, however, 
this term gives rise to spin relaxation in any real conduc- 
tor due to momentum scattering. Electron spins precess 
around an effective magnetic field that changes, as the 
momentum changes, at each scattering event. An en- 
semble of polarized spins, initially oriented in the same 
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FIG. 1: (a) An example of a trajectory segment in a 1/im-wide wire with A=10^m due to small angle scattering. Because the 
channel width is much less than the mean free path, the time dependence of the x- and y-components of the velocity, v^c and 
Vy respectively, are qualitatively different (panels (b) and (c)). (d) Power spectral densities show the square wave nature of 
ky oc Vy but no characteristic frequencies for oc v^. 



direction but following different random trajectories, will 
be distributed randomly around the Bloch sphere after a 
relaxation time, Tgr- This relaxation process is known as 
the Dyakonov-Perel' (DP) mechanism.'^' 

External magnetic fields have a strong effect on DP 
relaxation. These effects can be quite complicated when 
both orbital and spin effects are included. In this pa- 
per we discuss only the case of in-plane magnetic fields, 
which give rise to Zeeman splitting but not to Landau 
quantization or cyclotron motion. When Bso is added 
to an in-plane magnetic field B^xt, it is the total effec- 
tive field Btot = Bso + Bext that sets the spin precession 
axis and precession time for DP spin dynamics. The re- 
laxation time, Tsr-, due to the DP mechanism has been 
calculatecPSl for disordered 2D systems with momentum 
scattering time Tp, giving 
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Here Tp corresponds to a mean free path A — TpVp. 

The monotonic dependence Tg^iBext) described by 
Eq. (Sj) does not hold in confined geometries, such as 
the channels studied here, where the mean free path and 
spin-orbit length are on the order of or greater than the 
channel width. -'^ It is the goal of this work to study 
Tsr (Bext) numerically in these cases. 



III. MODEL 

Semi-classical Monte Carlo simulations of DP spin dy- 
namics in a 2DEG channel were performed by calculating 
momentum-dependent spin precession along an ensemble 
of randomly-generated classical trajectories, ri{t), analo- 
gous to the calculations described in Refs. I10I12I131 In- 
stantaneous velocity was determined from Vi(t) = dfi/dt, 
giving momentum Tikiit) = m*Vi{t) for effective mass m* . 
Throughout this paper the magnitude of the velocity was 
Wi\ = vp — lO'^m/s, cor resp onding to electron sheet den- 
sity w 5 X 10i°cm-2.[IE31 

Each spin Si was initialized to lie along the external 
field, Si{t = 0) |j Bext- The spins evolved in time by 
precessing around the trajectory-dependent Btotit), cal- 
culated using Eq. ^ and ki{t). The ensemble-averaged 
projection of the spin on the initial axis was then calcu- 
lated as a function of time, P{t) — {si{t) ■ 5^(0)),;, and 
fit to an exponential decay model, P{t) = Poe"*^^^'^, to 
extract the spin relaxation time Tsr- Although the DP 
mechanism can give rise to an oscillatory behavior for 
P{t) when the external magnetic field is zero (data not 
shown), the oscillatory component in P{t) disappeared 
above B^xt ~ for the device parameters studied here. 

The trajectories, fi{t), were confined to 1 ^m wide 
channels (Fig. 1(a)), refiecting the devices used in Refs.[5l 
1141 The range of spin-orbit parameters and mean free 
paths explored in this work represent a broader range 
than would commonly be encountered in a transport ex- 
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periment. The channels were assumed to be infinitely 
long, and each trajectory started from the middle of the 
wire with a random initial velocity direction. (It was 
confirmed that initial conditions had no effect on the cal- 
culated spin relaxation times after averaging over an en- 
semble of trajectories.) 

Disorder was taken into account primarily through 
small-angle scattering, although Fig. [5] compares the ef- 
fect of various scattering mechanisms and of soft vs. hard- 
wall confinement. With the exception of Fig. [5) scat- 
tering from channel walls was assumed to be specular. 
The semiclassical approximation used in this simulation- 
classical trajectories with coherent spin precession-is 
valid when both the orbital phase coherence time and 
the momentum scattering time are shorter than the spin 
relaxation time, and when electron trajectories can be 
assumed to be independent of spin direction. The lat- 
ter criterion implies that the Fermi energy is much larger 
than the spin-orbit energy -ffso- This is a valid approxi- 
mation in n-type GaAs quantum wells but not in p-type 
samples or narrow-gap semiconductors.'^ 



IV. RESULTS 

Figure [T] shows a short segment of a trajectory de- 
fined by small-angle scattering with mean free path of 
A = 10/im; this level of disorder that is experimentally 
accessible in high mobility electron gases. The qualita- 
tively different characteristics of the x- and y-components 
of momentum (Figs.jljb) and [ijc)) highlight the impor- 
tance of device geometry in the low-disorder regime stud- 
ied here. Electrons bounce off the channel walls many 
times before their momentum is randomized, because the 
mean free path is much larger than the channel width and 
scattering from the walls is specular. 

The trajectories in such a system are characterized by 
rapid, nearly-periodic changes in the sign of the momen- 
tum transverse to the channel, ky, while the magnitude of 
the longitudinal momentum, kx, changes only diffusively 
over a longer timescale. The square- wave character of 
ky{t) can be seen in Fig. mc), and in its power spec- 
trum Skyif) (Fig. ijd)). Notice that the ky frequency 
spectrum is strongly peaked despite the random angles 
of electron motion in a typical trajectory (Fig.[lja)). The 
peak frequency in Sky{f) reflects an average over the ran- 
dom distribution of trajectory angles, where the bounc- 
ing frequency for angle 9 is ^'^'^2w^^ ^'^^ channel width 
w. 

Relaxation of spins that are aligned initially along 
Bext results from fluctuating fields transverse to B^xt- 
In the DP mechanism, those transverse fields are the 
momentum-dependent effective fields arising from spin- 
orbit interaction. The first-order component of the effec- 
tive magnetic field due to ky is always in the x-direction 
(Eq. ([2])), independent of the relative strength of Rashba 
and Dresselhaus terms in the spin-orbit interaction. Sim- 
ilarly, the effective field due to kx is always in the y- 
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FIG. 2: (a) Spin relaxation times for x- (crosses) and y- (cir- 
cles) oriented magnetic fields, for trajectories with the same 
disorder parameters as shown in Fig. 1, a = 3meVA and 
/3 = 7 = (see text). Sohd and short dashed hnes are pro- 
portional to Skx.y~^ as described in the text. Long dashed 
line shows functional form expected in 2D. (b) Relax- 
ation times for different hnear spin-orbit coupling strengths, 
denoted as (a,/3). At high field the relaxation time Tsr{By) 
is determined by (a — /3), while at low fields Tar {By) also de- 
pends significantly on (a + jS). (c) Relaxation time at the first 
local minimum (rmm, By ~ 7T) and local maximum {rmax, 
By ^ 14. 5T), plotted against q (/3 = 7 = 0). Throughout the 
low-spin-orbit range (a < 6meVA), r,^ oc a (dashed Unes). 



direction. Hence spins in an external field along x relax 
due to fluctuations in the motion along x; spins in a field 
along y relax due to fluctuations in the motion along y. 

The qualitatively different power spectral densities 
of the two momentum components, Skxif) and Skyif) 
(Fig. 1(d)), give rise to qualitatively different relaxation 
behaviors for spins in x- and y-oriented fields respectively 
(Fig. 2(a)). Tsr{Bx) increases smoothly with Bx (spins 
initialized along x), matching the Tsr oc B^^^ behavior 
expected at high field in 2D disordered systems (Eq. ([s])) 
despite the confinement to a micron-wide channel in the 
simulation. The sharp periodic dips in Tsr{By) are the 
BSR features that are the subject of this paper: the short 
relaxation time at these dips is spin resonance due to the 
peak frequencies in Sy{f). This resonance occurs when 
peaks in Sy{f) occur at the Larmor frequency of the ex- 
ternal field, /l = gfisBext/h. 

Figure 2(a) also compares the inverse of the spec- 
tral densities of the two momentum components to the 
spin relaxation times, Tsr{Bext), extracted from the sim- 
ulations. Clearly, Sk~^{f) is directly proportional to 
Tgr{Bext) when Sk{f) is evaluated at the Larmor fre- 
quency, /i. This is reminiscent of the nuclear spin re- 
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FIG. 3: Efltect of tliird-order spin orbit term (quantified by 7) 
on BSR features. Eacli solid-line curve represents a different 
7 G {0, 50, 100, 500}eVA^, vertically offset as shown. Cir- 
cles represent 7 = SOeVA"^ (value expected for GaAs 2DEGs), 
and fall on top of the 7 = curve. All data is for a = 3meVA, 
/3 = 0, and A = lOfim. 



laxation time, Ti, for nuclear magnetic resonance, where 
it has been shown that Ti^^ — {gfj,B/^)^SB^{fL), with 
Sb± (/l) representing the spectral density of fluctuations 
in the transverse magnetic field, B±, at the Larmor fre- 
quency of the static NMR field. ^'"^ Fluctuations in Bgo 
are proportional to Sk{f ) by Eq. (2), and it is these fluc- 
tuations that lead to relaxation in the present case. As 
seen in Fig. 2(a), the approximation Tsr{B) oc Sk^^{f) 
becomes significantly worse when the mean free path is 
longer than lO^m, perhaps because the approximation of 
exponentially-correlated noise in the NMR result breaks 
down. 

Because Tsr{By) depends on the x-component of -Bso, 
it is controlled by (a — /?) and is nearly independent of 
{a + (3) (Eq. ([2|). The accuracy of this approximation 
can be tested in the simulation by varying a and {3 in- 
dependently. As seen in Fig. 2(b), curves with identical 
{a — p) but different {a + /3) fall on top of each other for 
By > 3T. (The stronger dependence on {a + (3) at low 
field comes about because the direction of Btot fluctuates 
significantly when Bf.xt ^ Bso-) The dependence of r^r 
on a is shown in Fig. 2(c), holding /3 = for all curves: 
when examined for particular values of magnetic field, 
the relation Ts 



expected for 2D (Eq. (|3 
carries over to the channel data (Fig. 2(c)). 

The discussion thus far has ignored the O(fc^) term 
in Eqs. (1) and (2). Averaged over the Fermi circle, the 
strength of this term is '^'ykp. Using values for 7 reported 

in the literature (9 — 34eVA^)pl and 2DEG parameters 
vp = I X lO^m/s and \a\ + |/3| « 3meVA reported in 
Ref. |SJ the third-order spin-orbit field, Bio' , is an order 
of magnitude smaller than the first-order field, isio' . For 
this reason, the simulations presented in most of this pa- 
per set 7 explicitly to zero for ease of calculation. 
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FIG. 4: (a) Relaxation times for different mean free paths 
in a 1/im-wide channel. Left- and right-hand panels show 
mean free paths less than and greater than 10/im respectively, 
with the A = lO/xm trace appearing in both panels. Signa- 
ture of ballistic spin resonance is clearly visible even when A 
is as short as the channel width l/im. Features continue to 
sharpen (including faster spin relaxation at resonance) up to 
A = 10/im. For longer mean free paths, relaxation time at 
resonance rises again, (b) Effect of mean free path on re- 
laxation time at the first local minimum (rmin, By ~ 7T) 
in a 1/im-wide channel. Arrows denote the crossover value 
Amin = 2nvFTso discusscd in text. 



For significantly larger values of vp or 7, on the other 



hand, ijio"* is of the same order or larger than Bso''. Be- 
cause of the more complicated symmetry of i?so\ its ef- 
fect on BSR is not monotonic in 7. Figure |3| explores the 
role of Bso by raising 7 while holding vp = 1 x lO^m/s, 
a = 3meVA, and j3 = 0. The BSR dips disappear around 
7 « 300eVA^ where sio^ « 2 x Bil\ but then revive for 
larger values of 7. 

Significant changes in spin relaxation were observed 
when the overall magnitude of disorder (set by A) was 
changed (Fig. |4]). When A was much smaller than the 
channel width, resonant dips were absent. In that case, 
the bouncing frequency ceases to be a relevant parameter, 
as electrons seldom make it across the channel without 
scattering, and the 2D limit of Eq. ([3| is approached. 
The dips become deeper as A is increased, but reach a 
minimum value around 10/xm before rising again for even 
longer mean free paths. 

In order to understand this non-monotonic depen- 
dence, we study the A-dependence of Tsr at the first res- 
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onant dip, around By = 7T (Fig. 4(b)). Starting from 
very short mean free paths, Tgr reaches a minimum at 
'^min — '^T^vpTso, then riscs again for very long mean free 
path. The length-scale, Xmm, corresponds to the dis- 
tance an electron would have to travel in order for the 
spin to rotate by 2n due to the spin-orbit effective field. 

This behavior can be explained at a qualitative level 
by considering spin relaxation in a reference frame that 
rotates at the Larmor frequency. Working in this frame 
effectively removes precession due the external field, and 
it removes flips in Bso that occur at frequency /l due 
to bouncing between the channel walls. In other words, 
spin relaxation in the ballistic channel at the BSR con- 
dition is approximately mapped onto spin relaxation in 
a disordered 2D system at zero external magnetic field. 
In 2D at zero external field, one expects r^r ^ tIoI'^p = 
t'^^vf/X (Eq. (3)) to decrease with increasing A in the 
motional narrowing regime, i.e. for fast momentum re- 
laxation, Tp < Tso- In the ballistic limit, Tp > t^o, on 
the other hand, one expects Tsr to increase with A as 
Tsr ^ Tp = X/v p bec ause spins precess coherently be- 
tween scattering! ^° l ^^ l It is the crossover from motional 
narrowing to ballistic regimes that gives rise to the non- 
monotonic behavior of Tsr in Fig. 4(b). 

Finally, we show that the particular type of disorder 
used to generate trajectories, and the type of scattering 
off channel walls, has only a small effect on the simulated 
spin relaxation curves. Figure [5] shows spin relaxation for 
three different types of disorder: 

1. small-angle scattering. The direction of motion 
changed from timestep to timestep by a small angle 
that was Gaussian-distributed around zero, with 
standard deviation calculated to give the desired 
mean free path. This is believed to be the domi- 
nant sc atteri ng mechanism in high-mobility GaAs 
2DEGS.I™ 

2. large-angle scattering was implemented as a prob- 
ability for complete randomization of momentum 
angle at each timestep. The probability was calcu- 
lated to give the desired mean free path. 

3. rough potential walls Upon reflection off channel 
walls, the angle of reflection was randomly dis- 
tributed around the angle of incidence with a 
spread of (j)spec- 4' spec — corresponds to specu- 
lar scattering from channel walls. This effect is be- 
lieved to be weak in electrostatically-defined GaAs 
2DEG nanostructures, as shown by clear transverse 
focusing signals even up to hish order, which re- 
quire many specular bounces.fiSl 

Each curve shown in Fig. 5(b) corresponds to disorder 
from only one of the three mechanisms. The mean free 
path is A = 10/xm in each case, confirmed by monitoring 



the autocorrelation of kx{t). As seen in the figure, the 
simulated spin relaxation time depends only slightly on 
the precise model of disorder, despite the importance of 
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FIG. 5: (a) Cartoon trajectories corresponding to the types 
of disorder and scattering described in the text, (b) Spin re- 
laxation time for B^xt || y, calculated for four disorder models, 
each giving effective X—10jj,m. Solid: small angle scattering. 
Dash-dot: large angle scattering. Short dash: rough walls 
(non-specular scattering), plus small angle scattering giving 
A=100/im for trajectories that would otherwise not encounter 
a wall. Dotted: small angle scattering and soft walls, mod- 
elled by a carrier density that decreases smoothly to zero 
within Xdep ~ 150nm of channel walls. Inset: match be- 
tween the power spectral density and spin relaxation time is 
much worse in the case of large angle scattering, compared to 
small angle scattering (see, e.g., Fig. 2(a)) and other disorder 
models (data not shown), {a = 3meVA, 13 — ■y — 0) 



ballistic transport to the resonant dips in Tsr{By). Figure 
5(b) also compares BSR for the case of simple reflections 
from hard-wall channel boundaries to the more realistic 
case of soft walls, with a 150nm depletion length as might 
be expected in nanostructures defined by electrostatic 
surface gates. Small angle scattering is implemented to 
give A = 10/im in both hard-and soft-wall simulations. 
The difference between the hard- and soft-wall data is 
nearly indistinguishable, except for a small shift in the 
field at which the resonance dips occur. 
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